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A THEORETICAL ANALYSIS OF THE EFFECTS OF FUEL MOTION ON AIRPLANE DYNAMICS 1 


By Albert A. Schy 


SUMMARY 

The general equations of motion for an airplane with a num- 
ber of spherical fuel tanks are derived. The motion of the fuel is 
approximated by the motion of solid pendulums. The same type 
of derivation and equations are shown to apply to any type of 
fuel tank where the motion of the fuel may be represented in 
terms of undamped harmonic oscillators. 

Motions are calculated for a present-day high-speed airplane 
and a free-flying airplane model with two spherical tanks in the 
symmetry plane. These calculations show that the normal air- 
plane motion may be considerably modified and that residual 
oscillations may result. The ratio of the natural fuel fre- 
quency to the natural airplane f requency is shown to be the most 
important parameter for determining the effect of the fuel motion 
on the airplane motion. The stabilizing effect of turbulence in 
the fuel is discussed , and it is suggested that the stabilizing efect 
of artificially induced turbulence be investigated experimentally . 

INTRODUCTION 

Small-amplitude lightly damped lateral oscillations are a 
troublesome characteristic of certain high-speed airplanes. 
Several possible explanations for these oscillations, which are 
adequate in specific cases, have been offered. For example, 
reference 1 shows that nonlinear aerodynamic derivatives 
could cause such oscillations, and it has been shown that 
atmospheric turbulence is another possible cause. It has 
also been suggested that a possible cause of such oscillations 
is the motion of fuel in the tanks. In some recently designed 
airplanes the mass of the fuel relative to the airplane mass is 
much larger than has been common in the past; therefore, 
the effects of fuel motion can be expected to be relatively 
more important. In fact, in several cases baffling the fuel 
tanks was found to have considerable effect on the general 
handling qualities of the airplane and sometimes actually 
eliminated the troublesome lightly damped lateral oscilla- 
tions which had been present. 

An experimental investigation of the effects of fuel motion 
on the lateral motion of a free-flying airplane model is de- 
scribed in reference 2. The results indicated that the effects 
of fuel motion were noticeable and caused the lateral motion 
of the model to be very erratic. 


The present analysis treats each fuel tank as a pendulum 
oscillating in two degrees of freedom and applies Lagrange’s 
equations of motion to obtain the interaction between these 
pendulums and the airplane. Thus, for small motions the 
fuels are treated as simple harmonic oscillators. The results 
are applied to obtain the general equations of motion of this 
system and, in particular, the lateral motion of an airplane 
with internal fuel tanks in the plane of symmetry of the air- 
plane. Since the general solution of the equations is ex- 
tremely complicated, an attempt is made to evaluate the 
results by carrying out numerical calculations for specific 
cases. This approach is shown to be adequate in yielding 
the most general effects of fuel motion. 

The discussion of the numerical application of the equa- 
tions of motion to specific cases is given in detail after the 
derivation of the equations of motion. This discussion of 
results is understandable quite independently of the deriva- 
tion of the equations of motion. 

SYMBOLS 


A, Y, Z 

airplane stability axes with origin deter- 
mined by equations (13); also compo- 
nents of applied forces along these axes 

L, M, N 

components of applied moments about 
X - , Y - , and Z-axes, respectively 

i, J k 

unit vectors along X- ) Y-, and Z-axes, 
respectively 

x, y , z 

components of translational displacement 
of airplane 

r 

vector translational velocity of airplane 
(ix +jy + kz = iU+v ) 

U 

magnitude of steady-state velocity 

U, V, w 

components of disturbance translational 
velocity of airplane 

V 

vector disturbance velocity of airplane 
( iu -fjv -f- kw) 

It 

vector position of a point in airplane 

(iR x -fjR y -\- kR z ) 

It/ 

vector position of center of gravity of 
fuel in a particular fuel tank 

V 

total vector velocity of a point in airplane 


1 Supersedes NACA TN 2280, “A Theoretical Analysis of the Effects of Fuel Motion on Airplane Dynamics” by Albert A. Schy, 1951. 
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sideslip angle ^tan- 1 —J 

infinitesimal rotations of airplane about 
X-, Y-, and Z-axes, respectively 
vector rotational velocity of airplane 

components of angular displacement from 
vertical of line joining fuel center of 
gravity to tank center, taken on mutu- 
ally perpendicular planes; f positive in 
direction of positive roll and 77 positive 
in direction of positive pitch 
distance from tank center to fuel center 
of gravity 

vertical displacement of fuel center of 
gravity from equilibrium position 
number of fuel tanks 
mass 

total mass of airplane and fuel 
total moments and product of inertia of 
airplane about X- } Y-, and Z-axes 
rigid-body moments and products of in- 
ertia about axes through center of gravity 
fuel moments of inertia about f- and 
r7-axes through tank center 


fuel moments of inertia about f- and 
77-axes through fuel center of gravity 
nondimensional radius of gyration in roll 

J 7 

m t b 2 

nondimensional radius of gyration in yaw 




nondimensional product-of-inertia param- 


kinetic energy 
potential energy 
period 

time for exponentially damped or increas- 
ing oscillation to halve or double am- 
plitude, respectively 
time 

nondimensional time parameter 
acceleration due to gravity 


flight-path angle with respect to horizontal 


P 

S 

b 

Vb f 


D 

C L 


Ci 


c 


71 


Cy 



Subscripts: 

J 


a 

0 


air density 
wing area 
wing span 

lateral nondimensional mass coefficient 
/ m t 

V 


P Sb ’ Mi, ' _ P SbJ 


™r \ 

-pSb) 

differentiation operator 

trim lift coefficient / — c - os ^ 

l y™ j 

rolling-moment coefficient 
/Rolling moment\ 


~ P U 2 Sb 


) 


yawing-moment coefficient 
/Yawing moment \ 


hpU*Sb 




-in ~ . , /Lateral force> 

lateral-lorce coefficient / - 


V 2 


>U 2 S 


particular fuel tank, or summation index 
over fuel tanks (/= 1 , 2 , ... k) 
airplane without fuel 
initial conditions at £=0 
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DERIVATION OF EQUATIONS OF MOTION 

ASSUMPTIONS FOR DERIVATION OF GENERAL EQUATIONS OF MOTION 

As a first approximation, only the effect of the motion of 
the fuel as a whole is considered; that is, only the funda- 
mental mode of the wave motion is considered, and this 
mode is approximated by rigid-body motion. The main 
effect of the internal wave motion is to introduce damping 
into the fuel oscillation. This damping is caused by the 
conversion of kinetic energy into heat through the turbu- 
lence caused by the splashing of the fuel. A strictly analytic 
consideration of such damping effects is extremely difficult ; 
on the other hand, the damping caused by the viscous tan- 
gential forces between the fuel and the tank is completely 
negligible (see reference 3). The analysis of the problem is 
therefore confined to the motion with no fuel damping and 
the effect of the damping is considered in the discussion of 
the results. 

In a spherical tank the fuel can oscillate approximately as 
a rigid body if no splashing is assumed for small oscillations. 
The motion ma}^ be pictured as the “rocking” of a spherical 
segment of constant shape. The restraining force of the 
tank, which always acts in a direction normal to the motion, 
is exactly analogous to the tension in a pendulum. Thus, 
the small motions of the fuel in a spherical tank may be rep- 
resented quite well by the well-known simple properties of 
small pendulum motions. This approach is used in the 
mathematical analysis of the problem. 

The effect of aspherical tank shape can be approximated 
by replacing the tank by an equivalent harmonic oscillator 
with an arbitrary amount of turbulence damping added even 
for small motions. For example, the representation of rec- 
tangular tanks as harmonic oscillators is discussed in refer- 
ence 4. Thus in this case also the most general effects of 
the fuel motion on the airplane motion should be qualita- 
tively obtainable by this type of analysis. 

The effects of large-amplitude fuel motions will be dis- 
cussed qualitatively after the discussion of the results of the 
mathematical analysis. As usual in stability analysis all 
motions are assumed small and second-order terms are 
ignored. 


DERIVATION OF GENERAL EQUATIONS OF MOTION 

With the preceding assumptions the physical problem can 
be considered as the interaction between two or more rigid 
bodies, namely the airplane and the several fuel pendulums, 
with each fuel pendulum considered as suspended from the 
tank center. The only potential energy considered in the 
system is that of the pendulums. If the inertial character- 
istics of the airplane and the fuel are known, the kinetic 
energy of the system can be obtained from the translational 
and rotational velocities of the airplane and the fuels. With 
this information the interactions in the system can be obtained 


by using Lagrange’s equations of motion in the form (see 
reference 5) 


dE k 


d /dEA 
dt\dqi ) dq t 


d E v 


c)(Z i 


Qi 


(*= 1 , 2 , . . . n) ( 1 ) 


where q t is one of the n generalized coordinates of the system 
corresponding to the n degrees of freedom, is the corre- 
sponding velocity, and Q t is the corresponding generalized 
force. The qi will be lengths and angles and the correspond- 
ing Qi will be forces and moments, respectively. 

The airplane itself introduces the customary six degrees 
of freedom, which are the three displacements of the airplane 
system along axes fixed in the airplane ( x , y , z) and the cor- 
responding angles of rotation of the airplane about these 
axes (<£, 0, i). For small displacements, the pendulum 
motion can be described by two angles f and tj since the 
vertical motion can be neglected (see fig. 1). The angle f is 
measured from a vertical line through the tank center to the 
projection of the line joining the tank center to the fuel 
center of gravity on the vertical plane parallel to the F-axis 
and tj is the corresponding angle in the vertical plane parallel 
to the X-axis. For small angles, f and 7? may be represented 



Vertical 


Figure 1 .— Decomposition of horizontal fuel motion in terms of angles f and r\ In vertical 

planes. 
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as in figure 1. In effect this figure makes use of the fact that 
small angles may be added vectorially. When the two addi- 
tional coordinates f and 77 are used to describe the pendulum 
motion, the whole system has two additional degrees of 
freedom for each fuel tank. 

Expressions must be obtained for E k and E p in terms of 
the coordinates of the system and their time derivatives in 
order to use equation (1). The only potential energy is that 
of the fuel pendulums, which can be written as follows for k 
fuel tanks: 

k 

E p = j: m fh f Q ( 2 ) 

/= 1 

For the height of the center of gravity in each fuel tank, as 
can he seen from figure 1, 

h = h{ + h v ~l — l COS f COS 7) — (l — \ f 2 )(l — ~ r)J 

or 

(3) 

Note that the vertical displacement h is of second order in 
the small quantities 77 and f. This fact justifies the previous 
statement that the vertical displacement could be neglected 
in describing the pendulum motion only by the two coor- 
dinates 77 and f. As might be expected, equations (2) and 
(3) indicate that each fuel pendulum is being considered as 
an undamped oscillator with two degrees of freedom in a 
horizontal plane. 

The kinetic energy of the total system can be written as 
the sum of the kinetic energy of the airplane and the kinetic 
energies of the fuels. Also the kinetic energy of each rigid 
body can be expressed as the sum of the translational energy 
of the mass moving with the velocity of its center of gravity 
and the rotational kinetic energy of the mass about its center 
of gravity. Thus, when the inertial characteristics of the air- 
plane and the fuels are known, the kinetic energy can be 
obtained as a function of the generalized coordinates and 
velocities if the translational velocity of each center of gravity 
and the angular velocities of the airplane and fuels about 
their respective centers of gravity can be expressed in terms 
of these generalized coordinates and velocities. 

In order to obtain the required expressions for these veloci- 
ties, a system of axes fixed in the airplane with the W-axis 
along the steady-state velocity at t=0 is used, as is customary 
in stability analysis. For the present the origin of the coor- 
dinates will not be specified. However, these stability axes 
are not inertial axes and Newton’s second law applies only in 
an inertial system of axes. The inertial axes may be taken 
as axes fixed in the earth. Then in the equations of motion 
the velocities and accelerations must be measured with re- 
spect to the earth, and their expressions in terms of com- 
ponents in the moving airplane axes may be obtained as 


shown in reference 6. These expressions will give the kinetic- 
reaction forces, which for the case of a rotating system are 
often referred to as “gyroscopic” forces. For the velocity, 
referred to the inertial system, of any point defined by the 
vector It in the airplane axes (in particular, for the centers 
of gravity previously discussed), 

V=r+coXlt-}-R==iU-\'V-{-(A)XR-\-R (4) 

where all vectors are given in terms of the airplane axes and 

r = ix+jy + kz = iU + 1 ? 

is the velocity of the origin of the airplane system with respect 
to the earth, while v and a > are the translational and rota- 
tional disturbance velocities of the airplane axes. 

Equation (4) may now be used to express the inertial 
velocities of the airplane center of gravity and the fuel 
centers of gravity in terms of the generalized coordinates by 
inserting for R the values R a and R/ , where R a is the vector 
position of the airplane center of gravity and R/ indicates 
the vector position of any particular fuel center of gravity. 
The vector R a is constant; therefore, R a =0. To obtain It/ % 
note that to first order 

R/~R f —il f sin (y o +0— 1 ?/) — jlfttf— 4> cos y 0 —\p sin y 0 ) + 

kl f cos (yo+0— *?/) 

R/ttR f —il f [sm y o +(0—rif) cos To]— .#/(?/— 4> cos y 0 — 

\p sin y 0 ) + A://[cos To— (0— V/) sin To] 

where R f is the fixed position of the tank center. Since To is 
constant, to first order 

R/~il f (rif—6) cos To jl/itf — 0cos To— t sin y 0 ) + 

sin To 

Again keeping only first-order terms leads to the following 
equation: 

u)XR/~o)XRf+ ilfi cos To — jl f (4 > cos To+ ^ sin to) + kl f 6 sin To 

Now combining the last two equations gives 

u)XR/+li/~uXRf+i'n f If cos To— jl/tf+kl f ri f sin y 0 (5) 

This equation shows, as could be predicted physically since 
no viscous force is assumed between the tank and the fuel, 
that the airplane rotation affects only the motion of the tank 
center. 

From equations (4) and (5) the necessaiy translational 
velocities can be obtained for the translational kinetic ener- 
gies. The rotational velocity of the airplane is simply w. 
The spinning motion of the fuel about the vertical axis is 
ignored; then, the rotational velocity of the fuel may be 
given by the components 77 and t • The two corresponding 
horizontal axes of rotation through the fuel center of gravity 
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are principal axes of the spherical segment of fuel; conse- 
quently, no product-of-inertia terms occur in the fuel rota- 
tional energy. Also, since the airplane center of gravity is 
in the airplane symmetry plane, Iyz=Ixy= 0 and only the 
Ixz product of inertia will appear in the airplane rotational 
energy. 

B y use of equation (4), the airplane velocity can be shown 
to be 

V a =i(U -]-U-\-dRz a — 'I'RyJ +j (0 +’/'^ar o ~0^2 a ) + 

k (w + <t> R Va — 6R Xa ) (6 a) 

When equation (5) is substituted into equation (4), the 
velocity of any particular fuel center of gravity is 

V f =i(U -\-u-\-6R Zf — iR Vf +vA c ° s To ) ~\~j (y-\-\!/R Xf — <t>Rz f — 

t/f) ^k(w-\-(t)Ry f —-dR Xf -\-rjflf sin 7 0 ) (6b) 

If V is the magnitude of the translational velocity of the 
center of gravity of a rigid body, /' is its moment-of-inertia 
tensor about the axes through its center of gravity, and co 
is the rotational velocity of the rigid body, the kinetic 
energy is 


when the components of the velocity (or any vector) are 
taken in the rotating system: 

-r—r+.Xr 

Thus the gyroscopic acceleration acting on the whole rotating 
system is 

c oXr =uX(iU+v)^(wXi)U=:jtU—kdU 


The effect of this acceleration can be brought into the equa- 
tions of motion by considering the inertial reaction of the 
total mass to this acceleration as an additional applied force. 

k 

If the total mass is m t = m a -\-^2m f , this reaction has the 
following components: 


Y' = — m t U\l/ 
Z' = m t Ud 


(9) 


In addition there is the inertial reaction torque M/ oil the 
fuel; this torque acts about the tank center and is caused by 
the acceleration of the tank center. For each tank, the 
vector reaction torque is 


E k =—m V 2 ^-^I x & x y' u v 2 -\-I z' — 2 1 yz u v u z 

2lxz'u x a) z 2 1 xy' ^ x ^v) (7) 

Thus, for the kinetic energy of the airplane, substitution of 
equation (6a) into equation (7) gives 

E k =^[(U+u+dR ! -iR v y+(v+tR x -*R: a y+ 

(w + 4>Ry a - 6R Xa y\ + I -f<i> 2 + I -^6 2 + I fi 2 -Ix Z '4>4' 

(8a) 

and, for the kinetic energy of each fuel, substitution of 
equation (6b) into equation (7) gives 

E kf = — [ ( 1 7 + u + 6R Zf — $Ry f + v f l f cos To ) 2 + 

(y + \l/R Xf — <t>R Zf — l flf ) 2 + ( W + <t>Ry f ~ 0R Xf -\- rjflf sin Yo) 2 ] + 


M/ = (R/ — R f ) X m f U(~j^ + kd) 

= m f U[(—il f sin t 0 + kl f cos y 0 ) X (— + kd) + 
Second-order terms] 

~ m f l f U{it cos to +ji sin y 0 +k\l/ sin 7 0 ) 

Since M^={M/) X cos 7 0 +(M /) 2 sin 7 0 and M/ = (M/) y , 

M/ = m f l f U \j/ (cos 2 7q + sin 2 7 0 ) = m f l f Ui ( 10 a) 

M/ = m f l f UQ sin y 0 ( 10 b) 

The forces and moments in equations (9) and ( 10 ) must 
be added to the weight and aerodynamic forces to obtain 
the Qi in equation ( 1 ). 

For convenience, the results of equations ( 2 ), (3), and ( 8 ) 
are as follows: 

Ep=~ gUm/lfitZ+vf) (1 la) 

Z f=i 


~2~Z' 2+ ~2 


V 


vr 


( 8 b) 


k 

For the total kinetic energy, E k =E k +^E k ; therefore, 

a /= 1 f 

equations (2), (3), and ( 8 ) may be used in equation ( 1 ) to 
obtain the ( 2 it + 6 ) equations of motion. However, it must 
again be recalled that the coordinate system is rotating. 
The whole system is therefore subject to an additional 
gyroscopic acceleration since the time-derivative operator 
contains an additional gyroscopic term (see reference 6 ) 


E k =^ [(U+u+6R ! -iR v ) 2 Hv+'i'R*-4>R* a ) i + 

(w+iR.-eRtf + 'jb [( U+V+9B,- 

\j/ R Vf -\-y)flf cos 7 o ) 2 +(^ ~\r'l'Rx f ~<l>Rz f ~ 

(w + d>Ry— 6R Xf + ilflf sin y 0 ) 2 ] |+^- 4> 2 + 

( 11 b) 
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Equations (9), (10), and (11) may now be used in equa- 
tion (1) to obtain the equations of motion. It should be 

d d d j d d ^ 

, and 7 ^= 5 - — For ex- 

ou oy ov oz ow 


. d d 
noted that = 5 - 777 - 1 — 5 

ox o(u +u) 


ample, to obtain the equation of motion in the z-direction, 
note that 

ZE k 


dx 


-^=0 

oz 


Then, equation (1) may be written as follows: 


<j_ / / ^>Ek \ 

dt\dxj dt\du ) 

= m a (u+dR 2 - iR Va ) + Xj m r(u + 0R Z — iR Vf + 
i )A cos 7o) 

= m,u + 6 (m a R Za + m/Rz ,) — ^( rn a R Va + 

S”iA / )+S’i/- m A cos 7 o=-X' (12) 

The position of the origin of airplane coordinates has not 

yet been specified. Equation (12) and the similar equations 
obtained in the other degrees of freedom suggest that the 
position of the origin be determined by the following three 
conditions: 

m a Rx a +J2 m f R X/ =0 (13a) 

m a R v +^2m f R v =0 (13b) 

a /= i 1 

m a R z +^m f R z = 0 (13c) 

° /= i 1 


Equations (13) imply that the origin is at the position of the 
total center of gravity when the fuel mass is treated as being 
concentrated at the tank center. This choice of the origin 
greatly simplifies expressions such as equation (12). The 
physical reason for this choice is again that the fuel does not 
rotate with the airplane; thus, a force acting on a line through 
this point, the center of gravity where the fuel reaction is 
assumed concentrated at the tank center, will produce no 
rotation of the airplane. 

The following substitutions will also greatly simplify the 
writing of the final equations of motion: 




/= 1 


- (14a) 


I z ^I/ + m4(R x y+(R v y}+^m f [(R Zf y+(R Vf y] 

k 

I = I xz +w a iL +]C m f R x R z 

a a f=\ f f 


Note that the quantities defined by equations (14a) are the 
total moments and product of inertia about the origin of the 


airplane coordinates when the fuel mass is assumed to be 
concentrated at the tank center. Finally, the necessary 
moments of inertia of each fuel pendulum about the tank 
center are 

h,=I r/+™A 2 




(14b) 


Without loss of generality then, equations (13) and (14) 
are used in the equations of motion obtained by substituting 
equations (9), (10), and (11) into equation (1). The general 
equations of motion can now be given as follows: 


+ S cos y 0 =X 

/= i 

k 

m t (w—Ud) + ^2 sin yo =2 

/= i 

I Y 6—(m a R Xa Ry a +^ m f R Xf R v ^ <£- 

^77 l a R Z( R Va -\- ^ j 'M'fRzfRv'f^\p ~ b 

^ (/? 2/ cos Yo Rx f sin Yo )m f lfij f =M 

f= i 

+ cos y 0 +{w—U6) sin Yo+ 

R Vl <j) sinYo+^ 2l cosYo— Rx x sin Yo) — 
Ry{P COS Yo] = 0 


- (15a) 


m 


t(v + 1 \p) — X! ^ 


/= 1 


I zX Ixz<t> (jft'aRy a Rz a ~\~ ^ j m f R Vf R z ^e 

k 

m fh{Rx f ^ f-\-Ry f V f COS Yo )=N 

I x <t> — Ixz'p — (^n a Rx a Ry a + Xj m fRx f R Vf ^ 0 + 

fc 

'y ' ^f^f{Rzf^f-\-RvfVf sin Yo) =jL 

i “ v — t ip+Rx l ii'—Rz l <t>) = 0 


(15b) 


In equations (15a) and (15b) only the fuel equations for the 
first fuel tank have been written. In each set there are k 
similar fuel equations. As has been previously stated, the 
forces on the right-hand sides of these equations are the 
applied forces and the weight and aerodynamic forces. 
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SIMPLIFYING ASSUMPTIONS 

The equations of motion have been separated into what 
would generally be considered the longitudinal motions, 
equations (15a), and the lateral motions, equations (15b). 
In the ordinary six-degree-of-freedom case, as can be shown 
from considerations of symmetry, no cross-coupling terms 
exist between these motions in the aerodynamic forces (see 
reference 7). Although such terms are known to exist in 
practice, they are small and generally neglected. However, 
many cross-coupling terms occur between equations (15a) 
and equations (15b) because of the fuel motions, even when 
the aerodynamic coupling forces are ignored. The magni- 
tudes of these fuel forces can be seen to depend on the masses 
of the fuels, the vector positions of the fuel tanks, the “ pendu- 
lum length” (i. e., the radius of the tank and the height of 
the fuel), and the accelerations involved in the fuel motion. 
Increasing the magnitude of any of these parameters will 
increase the effect of fuel motion. For most airplanes these 
fuel forces will be relatively small, but the present investiga- 
tion is primarily concerned with all first-order fuel effects. 
In any particular case, the actual magnitudes of these forces 
may of course be obtained by inserting the values of the 
previously mentioned parameters. 

Symmetrical fuel distribution. — Some of the terms in 
equations (15a)— for example, the <j> and \j/ terms in the 6 
equation — are essentially product-of-inertia terms arising 
from unsymmetrical distribution of the fuel about the y = 0 
plane (i. e., the plane of symmetry). That this is so can be 
seen if the fuel tanks are assumed to be distributed symmetri- 
cally with respect to the symmetry plane; that is 


S ™fRx,R Vf =ib rn f R z R y ==iti m f R v = 0 

/= i f= i f= i 


and from equations (13) 


ff „=0 


(16a) 


(16b) 


Therefore, these terms vanish for symmetrical fuel distribu- 
tions. In most cases the fuel will be symmetrically distrib- 
uted, and substitution of equations (16) into equations (15) 
yields the equations of motion for symmetrical fuel distri- 
bution: 

k > 

rn f lfijf cos y 0 =X 
/= i 


m t (w—U6) + ^2 Mflfirif sin y Q =Z 


/= i 


I Y 0 + lb cos to —Rif sin To 

Inji + mMgrjx+U COS To + (^— US) sin To+ 
R Vi (j> sin To+flCfr*, cos y 0 —R Xl sin To) — 
R v jp cos To] = 0 


►(17 a) 


m t (v+Uf )— S m Mr= Y 

/= i 


I z \js I xz<t> — ' 'ffl/lf(Rx r tf~\~ Ry r Vf COS To) — X 

/= l 1 


Ix<t> Ixzip- ^ jlfifi-z f~\~ Ry f Vf sin To) — L 


f= i 


f (17b) 


IiJi+rniliigti—v—Ut+Rx^—Rz^^O 


In equations (17) even though the terms arising from unsym- 
metrical fuel distribution have vanished, some cross-coupling 
terms still remain between equations (17a) and (17b). These 
terms occur in the 77 equations of set (17a) and in the </> and 
equations of set (17b). The significance of these terms is 
evident since each contains a factor R Vf . Thus, these terms 
arise when the airplane has fuel tanks with centers not in 
the plane of symmetry, even though they are symmetrically 
distributed with respect to this plane. For example, they 
would arise for wing- tip tanks. Physically, these terms 
clearly give the interaction between the longitudinal fuel 
motion r\ f in the wing-tip tanks and the airplane rotation 
about the vertical axis, which consists of the lateral motions 
</> and \p. For example, assume for simplicity that y () = 0 ; 
then, a yawing acceleration of the airplane will cause a lon- 
gitudinal fuel acceleration r\ f in the wing-tip tanks, and vice 
versa. 

From this discussion the 77 motion appears to couple the 
lateral and longitudinal airplane motions even for the per- 
fectly symmetrical fuel distribution described by equation 
(16a). However, the fact that this coupling does not occur 
can be seen by considering any pair of symmetrically placed 
and loaded tanks. Designate the 77 motion in this pair of 
tanks by 77 1 and rj 2 . Then, the 77 equations in equations (17a) 
show that a longitudinal horizontal acceleration of the tank 
gives rise as expected to 77 accelerations. Since the system 
is linear, this portion of the 77 motion may be considered 
independently, and because of the symmetry of the two 
tanks it is seen that vi=V 2 for the portion of the 77 motion 
arising from the longitudinal motion. Therefore, in the 
and <t> equations of set (17b) the effects of this 77 will vanish 
since R Vl =—R V2 . In a similar manner the laterally caused 
77 motion can be shown to have no effect on the longitudinal 
motion. Essentially the argument is that the 77 motion for 
each pair of tanks can be split up for perfectly symmetrical 
fuel distribution into symmetrical and antisymmet rical 
motions. The symmetrical portion of the 77 motion for each 
pair of tanks couples only with the longitudinal motion; the 
antisymmetrical 77 motion couples only with the lateral 
motion. Thus, for perfectly symmetrical fuel distributions 
the 77 equations of set (17a) could be combined with set (17b), 
only lateral degrees of freedom in the 77 equations being used 
(since the symmetrical portion of the 77 motion is of no 
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interest) ; or the rj equations may be used as shown with 
equations (17a), the lateral degrees of freedom in the ?/ 
equations being ignored. 

Fuel tanks centered in symmetry plane. — Many airplanes 
have large internal fuel tanks which are centered in the 
airplane symmetry plane. For such airplanes the equations 
of motion may always be separated into independent lateral 
and longitudinal modes since R Vf = 0. Using this value in 
equations (17) gives the equations of motion for tanks 
centered in the symmetry plane: 

k " 

V, 171 cos y 0 =X 
f= 1 

k 

mt(w—Ud) + y' l m f l f r\ f sin y 0 =Z 

/= i 


k 

IyO+'ZZ (Fz, c°s To Rx f sin To )mjl f rj f =M 
/= i f 1 

Ivji + mMgVi+ucos To + (w—U6) sin y 0 + 

( R Zl cos To —R X{ sin 7 0 )£] = 0 


r (18a) 


m t (v + U\p) — ^2 m f l f l f =Y 

/= i 


I zip I XZ<t> TUflfRx’i f — N 


/= 1 


lx <t> Ixz^~\- IzZ TUflfRz if — L 

f= 1 f 


► (18b) 


Ititi+militgti— v — U\p—R Xl ^—R Zl <j>) — 0 


Each set of equations, (18a) and (18b), contains (£+3) 
variables; the two sets can be seen to be independent of each 
other since the f motion of each tank couples with only the 
lateral motion and the 77 motion couples with only the longi- 
tudinal motion. 

For the case of a single fuel tank at the airplane center of 
gravity the modification of equations (18) is obvious. Then 
R a =zR f =0 } and all coupling between the rotational motion 
and the fuel motion vanishes; that is, all the fuel terms in the 
rotational equations vanish and all rotational terms in the 
fuel equations vanish. For an aspherical tank the rotational 
coupling in this case will be small. 

LIMITATIONS INHERENT IN THE APPROXIMATIONS 

Before proceeding to the application of equations (18b) 
it is appropriate to consider somewhat more explicitly the 


assumptions involved in the indiscriminate dropping of all 
second-order terms which appeared during the derivation 
of the equations of motion. In this connection the correc- 
tion, arising from the airplane accelerations, to the constant 
acceleration field g involved in the pendulum potential en- 
ergy should be considered. The assumption which is implied 
in neglecting these accelerations is that the accelerations of 
the tank centers are small with respect to g. 

If this and previous approximations are considered, it 
can be seen that three essential assumptions were made in 
dropping second-order terms: 

(1) The fuel and airplane angular displacement variables 
are small enough so that the angle approximates its sine. 
However, this approximation sometimes took the form that 
the angle was much less than 1 radian. 

(2) The disturbance velocities are much less than U> and 
products of the linear or angular velocities can be ignored. 

(3) The accelerations of the tank centers must be small 
compared with g. 

Strictly speaking then, the statement that the equations 
of motion (15), and also the simplified equations, are ac- 
curate equations of motion to first order is to be taken to 
mean that the motions to which these equations apply are 
restricted by the preceding three conditions. Thus, the 
equations would appear to remain accurate at least at the 
beginning of a disturbance. Moreover, when the motion 
becomes large enough that these assumptions break down, 
the fundamental physical assumption that the fuel may be 
considered to move as a rigid body also breaks down; there- 
fore, nothing can essentially be gained by keeping higher- 
order terms in the mathematical expressions. 

Since the pendulum motion is little changed even up to 
angles of 30° to 40°, it could be expected that aside from 
splashing effects these equations should remain a good to 
fair approximation even at such angles. On the other hand, 
even the splashing effects, although they would introduce 
some damping and change the inertial characteristics of the 
pendulum somewhat, could certainly not be expected to 
cause the general assumptions to break down completely 
for fuel motions up to angles of 30° to 40°. Therefore, the 
equations of motion derived are assumed to present a fair 
picture of the disturbance motion even up to fuel displace- 
ments of this magnitude. 

APPLICATIONS TO SEVERAL CASES 

NON I) I MENSION AL EQUATIONS FOR TANKS IN SYMMETRY PLANE 

The equations of motion (18b) have been applied to the 
lateral motion in several cases with two fuel tanks in the 
plane of symmetry. In these cases the lateral motion can 
be considered independently. The applied forces are the 
weight, the usual aerodynamic forces linear in the disturb- 
ance velocities, and any disturbing forces that may be 
present. In order to put the equations in nondimensional 
form, the nondimensional lateral airplane equations are used 
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as obtained in reference 8. The fuel equations are made 
nondimensional by making the standard transformation to 
nondimension al time derivative, as in the airplane equations, 

7YI l ^ 

and then by dividing through by • The resulting 

nondimensional expressions in the following equations are 
defined in equations (20): 


(yUbl ) — ~ C Y ^ p +(^m,D — i C L tan Yn) i— 2 C L <t>— D 2 { 




Cy 


(19a) 


~\ C„ r D ) i+(^»K xz D*-\ C n D ) <t>- 

(19b) 

-lc,'p+(» t K xz D>-±C, r Dy+(^K x W*-±C lp Dy+ 

Zj i Z 2 -Z? 2o 1 

m», -yZ) 2 f,+ MS2 -^D^ 2 =^c, (19c) 

£> 2 +p)^+^ 21 D 2 «+(^, 2 D 2 +^)fi=0 (1 9d) 

-Dp-(jf- D*+d) Z? 2 *+(Wz? 2 +G^ 2 =0 (1 9e) 


where 




( 20 ) 


The derivatives C Y and C Y were assumed to be zero. 

■* p * r 


METHODS OF SOLUTION 

In the present case the two fuel degrees of freedom intro- 
duce two additional oscillatory modes into the characteristic 
solution, in addition to modifying the original airplane mode. 
The motion will therefore be a combination of three oscilla- 
tions (aside from the less important exponential modes), but 
just knowing the three oscillatory roots is insufficient to indi- 
cate the type of motion since the relative magnitudes of the 
oscillatory modes must also be known. For this reason 
motions must be calculated in order to see the actual effects 
of the fuel motion. However, in several cases the charac- 
teristic roots were also found in order to facilitate the inter- 
pretation of the motions. These cases will be discussed 
subsequently. 


The most convenient method for the analytical solution 
of a set of linear ordinary differential equations such as 
equations (19) is probably the Laplace transform method 
(see reference 8). However this method is extremely cum- 
bersome and difficult to check since it involves the expansion 
of fifth-order determinants in which the elements are often 
quadratic functions of the characteristic root. Therefore, it 
seemed preferable to use some step-by-step method which 
would be more amenable to machine computation. 

Reference 9 gives a matrix method for getting the step-by- 
step solution of a set of linear ordinary differential equations. 
When applied in the present case to equations (19) this 
method results in a simultaneous solution for the motion in 
each of the five degrees of freedom and also for the motion 
in D(j), D\p, D£i, and D{ 2 - The calculations were carried out 
on the Bell Telephone Laboratories X-66744 relay computer 
in use at the Langley Laboratory. The essential details of 
the method are given in appendix A. 

SOLUTIONS FOR SEVERAL CASES 

The two basic cases for which motions were calculated 
were case A, a present-day high-speed airplane with two fuel 
tanks satisfying the conditions for equations (19), and case 
B, which corresponds essentially to case B of the model used 
in reference 2. The essential parameters for these two cases 
are given in tables I and II. Table II gives the conditions 
for case A when both tanks are one-half full (Ai) and when 
the fuel height equals one-half the radius (A 2 ) and for case 
B when the fuel heights in both tanks are 2 inches (B 2 ), 3 
inches (B 3 ), and 4 inches (B 4 ). 

In case A the tanks are spherical, somewhat over 4 feet in 
diameter, and centered on the body axis approximately 4 
feet in front of and behind the airplane center of gravity. 
The flight conditions are given in table I. The fuel weight- 
in the half -full condition is approximately 25 percent of the 
total weight. 

In case B the tanks are spherical, centered in the plane of 
symmetry slightly less than 5 inches below the airplane axis 


TABLE I 

STABILITY" DERIVATIVES AND MASS AND GEOMETRICAL 
CHARACTERISTICS FOR TWO CASES CONSIDERED 


Parameters 

Case A 

Case B 

Cyo - 

-1.042 

-0. 80 

w 0 

C n « 

0.17 

•0. 17 

Cu . 

-0. 126 

-0. 14 

‘0 

Cn - ..... 

-0. 01552 

b -0. 040 

K ' n P 

Cl 

-0. 342 

-0. 30 

Cn - - 

-0.28 

-0. 16 

Cl . 

0. 0796 

0. 30 

Weight of airplane alone, lb — 

6,970 

11.25 

S, sq ft . 

130 

2. 67 

6, ft 

28 

4.0 

*y 0 , deg - - 

0 

-11 

p, slugs/cu ft 

0. 00136 

0. 002378 

U, ft/sec 

704 

61. 5 to 74. 25 


» In case Bu, C n ^=0.29. 

b Actually, slightly different values of C np were used for each of the subcases of case B. 
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and 4 inches in front of and behind the model center of 
gravity. The diameters are 8 inches and the total fuel 
weight in the half-full condition is approximately 46 percent 
of the total weight. The flight conditions are the same as 
in reference 2. 

Motions were calculated for certain subcases of the basic 
cases which were obtained by varying certain significant 
parameters. By comparing the resulting motions an attempt 
was made to evaluate the effect of varying such factors as 
amount of fuel, position of tanks, and relative natural fuel 
and airplane frequencies on the disturbance of the airplane 
motion caused by fuel motion. Also various initial condi- 
tions were considered to show the effect of initial conditions 


TABLE II 

FUEL DEPENDENT PARAMETERS 


Parameters 

A, 

A a 

P 2 

b 3 

B< 

Radius of forward 
tank, ft 

2. 15 

2. 15 

0. 333 

0. 333 

0. 333 

Radius of rear tank, ft 

2.12 

2. 12 

0. 333 

0. 333 

0.333 

Fuel height in for- 
ward tank, ft 

2. 15 

1.075 

0. 1667 

0. 250 

0. 333 

Fuel height in rear 
tank, ft 

2. 12 

1.06 

0. 1667 

0. 250 

0. 333 

R,,, ft.. 

3.5 

3.5 

0. 333 

0. 333 

0. 333 

Rt V ft 

-4. 1 

-4. 1 

-0. 333 

-0. 333 

-0. 333 

Rz v ft 

0. 0123 

0.0123 

0. 407 

0. 407 

0. 407 

R, 2 , ft 

-0. 0144 

-0. 0144 

0. 407 

0. 407 

0. 407 

Weight of forward 
fuel, lb 

1,480 

462 

1.51 

3. 07 

4. 85 

Weight of rear fuel, lb 

1,068 

334 

1.51 

3. 07 

4.85 

Total weight of air- 
plane and fuel, lb. 

9. 518 

7. 766 

14. 27 

17. 39 

20.95 

Cl — . 

n 19 

0. 19 

1.062 

1. 178 

1. 174 

Ix, slug-ft 2 ... 

1,360 

1,360 

0. 1081 

0. 1243 

0. 1428 

//. slug- ft-' 

7, 340 

7,708 

0 2()9S 

0. 2206 

0.233 

1X2, Slug-ft 2 

262 

277 

0 

0 

0 

ll, ft 

ii. 806 

1.45 

0. 225 

0.174 

0. 125 

h, ft 

0. 795 

1.43 

0. 226 

0. 174 

0. 125 

Jj, slug-ft 2 

85.2 

40.0 

0. 00315 

0. 00507 

0. 00668 

h, slug-ft* 

59.7 

28. 1 

0. 00315 

0. 00507 

0. 006G8 


on the resulting motion. In some cases an initial disturbance 
in sideslip was assumed, and in other cases an initial fuel 
disturbance was assumed. An initial sideslip of 5° and a 
fuel displacement of 10° were arbitrarily chosen as standard. 
Since the equations are linear, multiplying the initial dis- 
placements by a common factor simply multiplies the result- 
ing motions by the same factor. In this connection it must 
be emphasized that, if at any time the calculated motion in 
any degree of freedom becomes too large to satisfy the 
limitations previously discussed, the following motion is 
meaningless. For example, if a 5° displacement in P gives 
rise to a fuel motion much greater than 30° to 40°, the accom- 
panying P motion is meaningless because the assumption of 
small displacements is violated. However, if multiplying the 
fuel motion by some arbitrary factor, for example 2/5, will 
bring its peaks down to less than 30° to 40°, then the p motion 
resulting from an initial P disturbance of 2° can be obtained 
by simply multiplying the previous P motion by 2/5 also. 
The effect of large fuel displacements must be discussed 
qualitatively. 

The motion in sideslip and the motion of the two fuel 
pendulums in the various subcases are shown in figures 2 to 
12. Comments on these motions are presented to facilitate 
interpretation of the figures. The period of each fuel pendu- 
lum is called the natural fuel period. The period and damp- 
ing of the airplane, the fuel being disregarded, are called the 
natural airplane period and damping. 

Case A lt — The natural fuel periods for case Ax (half-full 
tank) are approximately 1.66 seconds and the natural period 
of the airplane alone is 1.40 seconds. Damping to half- 
amplitude occurs in 2 cycles. 

The motion in figure 2 is for initial p 0 — 2°. The early p 
motion seems to have more damping than the natural airplane 



Figure 2. — Present-day high-speed airplane with tanks one-half full (case Ai). Sideslip and fuel motions following initial sideslip, fJo—2°. 
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mode. The disturbance arising from the fuel modes is evi- 
dent after 2 cycles. That the irregular residual oscillation 
of amplitude to %° is essentially due to the fuel modes 
is seen from the fact that the dominant period in the later 
motion is approximately 1.6 seconds. Notice that in this 
case a 5° initial would almost immediately cause fuel dis- 
placements of over 80°, so that the following motion would 
be radically changed. 

The motions shown in figure 3 for (fi)o= (£ 2 ) 0 = 10° are 
quite regular and indicate one dominant mode in each motion. 
The fuel period is 1.6 seconds. The airplane period starts 
at 1.4 seconds, builds up to 1.7 seconds, and averages 1.6 
seconds. The amplitude of the sideslip motion is very small. 
The largest such motion which coidd occur for this type of 


disturbance would be for initial (fi)o= (£* 2 ) 0 = 30 ° and would 
give ft amplitude slightly more than 0.1°. 

The small amplitude of the sideslip motion in figure 3 was 
surprising. It was conjectured that for this fuel configura- 
tion the fuel displacement (fi) 0 =10° and (£ 2 ) 0 = — 10°, 
corresponding essentially to an initial yawing moment, might 
be more effective in inducing an airplane oscillation. (See 
fig. 4.) Apparently, this configuration is more effective 
inasmuch as the sideslip motion now builds up to an ampli- 
tude of approximately 0.4°. The energy necessary to induce 
this considerable “snaking” type of oscillation seems to be 
obtained initially from the rear-tank motion, which is in the 
proper phase relation with the sideslip motion to feed energy 
into it at the start of the motion. 




Fioure 3.— Present-day high-speed airplane with tanks one-half full (case A»). Sideslip and fuel motions following initial fuel disturbance, (i"i)3=(fOo=10°. 
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The fuel periods in the regular motion arc slightly over 
1.6 seconds. The airplane period increases from 1.5 to 1.8 
seconds and has an average period of 1.6 seconds. 

Case A 2 .* — The natural fuel periods in case A 2 are approxi- 
mately 1.52 seconds and the natural airplane period is 1.49 
seconds. Damping to half-amplitude occurs in 1.3 cycles. 

The motion in figure 5 is for initial /3 0 =0.5°. The early 
sideslip motion seems to be of greater damping than the 
natural airplane mode. The residual airplane motion arising 
from the fuel modes sets in very quickly and is a regular 
unstable motion of very large relative amplitude with a 


period of approximately 1.7 seconds. Both fuels start with 
a period of approximately 1.5 seconds, which increases to 
1.7 seconds. 

In figure 6, (f!) 0 = — 10° and (f 2 ) 0 = 10°. The sideslip 
builds up to a fairly regular oscillation of 0.4° amplitude 
with the period increasing from 1.5 seconds to 1.7 seconds. 
The fuel motion has a period somewhat under 1.7 seconds 
with amplitude quickly building up to the limits where 
splashing must become important. 

Case B 2 . — The sideslip motion shown for case B 2 appears 
to be a normal damped oscillation for the first 4 seconds 



Figure 4.— Present-day high-speed airplane with tanks one-half full (case Ai). Sideslip and fuel motions following initial fuel disturbance, (fi)o— 10°, (r 2 )o— — 10°. 
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(fig. 7, p 0 =5°), but then the peaks show a slight irregularity 
instead of damping smoothly. The period of the early side- 
slip motion appears to be somewhat over 0.9 second and the 
motion damps to half-amplitude in less than 2 cycles. This 
motion is very close to the undisturbed airplane model 
(period of 0.92 sec and damping constant of 1% cycles). 
The natural fuel periods are 0.61 second. The fuel motion 
is very irregular and obviously contains considerable amounts 
of at least two characteristic modes. 

Cases B 3 and B 4 . — The motion in cases B 3 and B 4 shown in 
figures 8 and 9, respectively, for 0o=5 o is veiy much the 
same as in the previous one, except that the disturbance of 
the airplane mode in the sideslip motion appears somewhat 
more pronounced as the amount of fuel increases. 

Figure 10 (case B 4 , (fi)o=(f2)o=10°) shows that the motion 
in sideslip resulting from the fuel displacement is much more 
irregular than in case A. The dominant mode corresponds 
to a fuel frequency, but apparently the airplane mode is 
present with considerable amplitude. The maximum oscilla- 
tions are approximately ±#°. The sideslip motion in this 
case was much more irregular than for the corresponding 
initial conditions in case A. It was conjectured that this 
might be caused by the fact that in this model both tanks 
are below the X-axis, so that the coupling of the fuel motion 
with yawing and rolling motions does not have the same 
phase relationship as in case A where one tank is above and 


one is below the X-axis. Therefore in case B 4a , shown in 
figure 11, the front fuel tank was assumed to be above the 
X-axis, all other conditions remaining as in case B 4 . In this 
case the general type of motion does seem to resemble that 
in figure 3. The sideslip, which builds up to ±M°, shows 
a snaking at the fuel frequency. 

In case A the fuel natural periods are very close to the 
airplane period. In case B, however, the fuel period is ap- 
proximately two-thirds of the airplane period. In case B 4& , 
shown in figure 12 for /3 0 =5°, the value of C np of the model 
has been arbitrarily changed to give the model a period very 
close to the fuel period of 0.66 second. Comparison with 
figures 2 and 5 shows that the motion in this case is very 
much like the motion in case A. 

Transverse accelerations. -In evaluating pilots’ reactions 
to snaking oscillations, the magnitude of the transverse 
accelerations involved in the oscillations has been found to 
be an important factor. Acceleration amplitudes above 
0.025*7 are found to be bothersome, and amplitudes above 
0.08*/ are considered very unsatisfactory. Calculations of 
the transverse accelerations involved in several of the previ- 
ous motions were carried out. The magnitudes of the accel- 
eration peaks in the residual oscillations were found to be 
approximately 0.04*7 to 0.05*7. The actual motions are not 
shown since all the airplane oscillations are essentially of the 
same type as the 0 motions. 



Figure 5.— Present-day high-speed airplane with fuel heights one-half the tank radius (case Aj). Sideslip and fuel motions following initial sideslip, 0o=O.5°. 
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Figure 7.— Free-flying airplane model with fuel heights one-half the tank radius (case B 2 ). Sideslip and fuel motions following initial sideslip, 0o=5°. 


DISCUSSION OF RESULTS 

The motions described in the preceding section are suffi- 
cient to give a fair picture of the types of possible fuel effects. 
Moreover, since each motion is just a superposition of the 
characteristic modes of the total system, these motions are 
often easier to understand if the characteristic roots are 
known. Physically it is clear that the characteristic modes 
will not differ much from the natural (uncoupled) modes 
when the interaction between the airplane and fuel is small. 
Comparison of figures 2, 5, and 12 with figures 7, 8, and 9 
indicates that the interaction between airplane and fuel is 
strongest when the frequency of the airplane is close to that 
of the fuel, as might be expected from comparison with the 
resonance phenomena exhibited by an oscillator driving a 
system at its natural frequency. For this reason the charac- 
teristic modes of the total system were calculated in cases 
A 2 and B 46 , where the frequency ratio between airplane and 
fuel natural frequencies was practically unity. The natural 
modes are given for purpose of comparison. The results 


are given as follows in terms of periods and times to halve 
or double amplitudes, in seconds: 

Case A 2 , natural modes: 

P a =lA9 T H = 1.91 

Pi=P 2 =1.52 

Case A 2 , total system: 

P a =1.47 T »= 61 

Pi = 1.29 jP^=1.34 

P 2 =1.67 T 2 = 4.56 

Case B 46 , natural modes: 

P a =0.66 P H =1.29 

P X =P 2 = 0.66 
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Figure 8— Free-flying airplane model with fuel heights three-fourths the tank radius (case Bj). Sideslip and fuel motions following initial sideslip, 0o=5°. 


Case B 46 , total system: 

P a =0.65 T v = 30 

P, = 0.54 T 2 = 1.5 

P 2 = 0.71 T h =0A8 

In these eases of large interaetion it is diffieult to identify 
one of the characteristic modes as the airplane mode. The 
characteristic mode in which the period is changed least 
from the natural airplane period has been called the airplane 
characteristic mode. However, in figures 5 and 12 this 
mode is not obviously the dominant one, as the airplane 
mode is in figures 7, 8, and 9. 

The most important effects to be noted in these particular 
cases, where the fuel and airplane frequencies are equal, are 
that a characteristic mode which is very lightly damped with 


a frequency close to the natural airplane frequency exists 
and that an unstable mode appears. In connection with the 
first of these effects, it would seem that, theoretically, 
certain initial conditions might be found that would excite 
mainly this lightly damped mode in the characteristic solu- 
tion for p, so that the resulting motion would be a typical 
snaking. Of course the required initial conditions might or 
might not be practical ones. 

The total characteristic modes were also calculated for 
cases B 4 and B 4a to investigate the changes in motion caused 
by a hypothetical shift of one of the fuel tanks. The natural 
modes and the characteristic modes of the total system are 
given as follows for comparison: 

Cases B 4 and B 4a , natural modes: 

P a =0.84 T h =1M 

P 4 = — P 2 — — 0.66 
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Case B 4 , total system: 

P a =0.88 7^=1. 91 

P^O.63 P h =5.45 

P 2 =0.53 ^=1.61 

Case B 4a , total system: 

P fl =0.86 7^= 1.14 

P x = 0.61 Ty= 19.5 

P 2 =0.60 P 2 =10.2 


It is interesting to note that both fuels in figure 10 and 
also in figure 11 seem to follow the more stable fuel mode in 
the part of the motion shown. The 0 motion in figure 11 
seems to show the effect of the unstable mode. It appears 
in this case that, when the fuel tanks are in front of and 
behind the center of gravity, the configuration with one tank 
above and one below the X-axis gives rise to an unstable 
mode; whereas the configuration with both tanks below the 
X-axis makes both fuel modes stable. 

A comparison of figures 2 to 4 or figures 5 and 6 clearly 
shows that the initial conditions can have a very important 
effect, since the least stable mode does not necessarily become 
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Figure 10.— Free-flying airplane model with tanks one-half full (case B<). Sideslip and fuel motions following initial fuel disturbance, (O)o=(f2)o=10°. 



Figure 11 -Free-flying airplane model with tanks one-half full and front tank assumed to be shifted above the X-axis (case B<.). Sideslip and fuel motions following initial fuel disturbance, 

(fi)o— (fa)o=10°. 


dominant for a long time. The 13 disturbance was chosen as 
a typical airplane disturbance. On the other hand the dis- 
turbance of the fuels as an initial condition would seem to 
be completely artificial. However, these motions are be- 
lieved to give a rough idea of the residual oscillations caused 
by fuel motion, at least insofar as magnitude is concerned, 
since, if the fuels were still displaced after the airplane motion 


had practically died out, the remaining motion might be 
considered to be the type caused by a fuel displacement. 
From this point of view figures 4, 6, 10, and 11 seem to 
indicate that residual oscillations of the order of magnitude 
of /° to y 2 ° might be expected in these cases. Actually, 
figures 10 and 11 would show oscillations of the order of 
%° for 30° fuel displacements. 
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Figure 12.— Free-flying airplane model with tanks one-half full and C n p adjusted to make natural fuel and airplane frequencies equal (case Ba). Sideslip and fuel motions following 

initial sideslip, 0o=5°. 


Because of the lengtliiness of the calculations, only cases 
Ai, A 2 , and B 46 were carried out far enough to show the 
residual oscillations following a 0 disturbance. It is evident 
that the motion in case B 4& (fig. 12) resembles the motion in 
case A (figs. 2 and 5) much more than it resembles the un- 
modified case B motion (figs. 7 to 9). The reason for the 
smaller relative fuel motion in case B 4& is probably the fact 
that the relative fuel mass is considerably larger than in 
case A. In case A 2 (fig. 5) the residual oscillation dominates 
the motion almost immediately. These results show that the 
importance of the residual oscillation depends mainly on the 
closeness of the natural airplane and fuel frequencies, that is, 
on the parameter which might be called the frequency ratio. 
Inasmuch as the previous discussion of the characteristic 
modes indicated that the frequency ratio was also the most 
important factor affecting the characteristic modes of the 
system, the frequency ratio generally can be seen to be the 
most important factor determining the disturbance of the 
normal airplane motion caused by the fuel. Moreover, case 
A indicates that for spherical tanks the fuel frequency may 


easily be of the same order of magnitude as the airplane 
frequency. Reference 4 indicates that the same is true for 
rectangular tanks and for arbitrarily shaped tanks of reason- 
able dimensions. Thus, even though the residual oscillations 
might occur at fuel frequencies, these frequencies would not 
be distinguishable from the normal airplane frequency in the 
cases where the fuel effect is most pronounced, since in these 
cases the frequency ratio approaches unity. 

The effects of unstable modes cannot be understood with- 
out considering the nonlinear effects due to splashing of the 
fuel. For linear systems the presence of an unstable mode 
would imply that the total system is unstable. It has been 
shown, however, that in an actual motion if the coefficient 
of the unstable mode in the solution for the airplane motion 
is very small compared with the coefficient of one of the 
stable modes, then the unstable mode will not appear in the 
early part of the motion. Now even in cases where the inter- 
action is weak, one of the fuel modes (with no natural damp- 
ing assumed) may be unstable. In such cases the unstable 
mode in the airplane motions will be relatively very small, 
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while the unstable mode may be dominant in one of the fuel 
motions. Then this fuel motion may become very large 
before the effect on the airplane motion can be observed. 
At this point splashing will set in and the equations of motion 
no longer hold. Actually the energy lost in splashing may 
damp out the fuel motion so that the effect of this fuel mo- 
tion, and therefore of the unstable mode, may never ap- 
preciably appear in the airplane motion. 

In case A 2 , however, the interaction is large, and figure 5 
shows that the unstable mode does soon become very important 
in the /3 motion. Here again the effects of large-amplitude 
fuel motion and, in particular, the nonlinear effects of fuel 
splashing should be considered qualitatively. In cases Ai 
and A 2 , because of the strong interaction, a sideslip of 5° 
woidd very quickly give rise to fuel displacements of the 
order of magnitude of 90° or greater. Physically, it is clear 
that the airplane motion is feeding energy into the fuel 
motion in this early part of the motion. Because of the 
large amplitudes of the fuel motion, a large component of 
this motion is in the vertical direction, so that a considerable 
part of the energy in the fuel motion will go into creating 
longitudinal airplane motion. Because of the symmetry con- 
ditions previously cited this energy will not be fed back into 
lateral motion, so that the fuel tends to stabilize the lateral 
motion by feeding some of the energy from the lateral into 
the longitudinal motion. Also the turbulence due to splash- 
ing will absorb energy which will then be lost altogether 
from the motion. 

One conclusion which can be drawn from this discussion 
is that the motion following large disturbances may be more 
stable than that following small disturbances. For example, 
the initial disturbances in figures 2 and 5 have been adjusted 
so that during the motion shown there is little splashing, and 
considerable residual oscillations are shown. But if the ini- 
tial disturbances in these cases had been 5° or more in 0, 
the energy lost in splashing in the early motion would pos- 
sibly be so great that the residual motion in /3 would be 
smaller than that shown here. Tliis might explain why some 
airplanes which definitely showed troublesome fuel oscilla- 
tions were reported to be more stable in conditions of large 
atmospheric turbulence than in slightly turbulent atmos- 
pheric conditions. 

A more important conclusion is that the fuel can be used 
to stabilize the airplane motion by introducing turbulence — 
for instance, by use of appropriate baffles. This can be seen 
by noticing that the early part of the 0 motion in figures 2 
and 5 is very stable. But if most of the energy fed into the 
fuel in this part of the motion were converted into heat 
through turbulence, then as has been pointed out this energy 
could not be fed back into the airplane motion and the 


residual oscillations would not appear. Since the amount 
of energy lost in turbulence cannot be calculated analytically, 
it would seem that an experimental investigation of the 
effects of honeycomb or other turbulence-inducing baffles 
on the airplane stability w^ould be desirable, especially in 
cases where the airplane and fuel natural frequencies are 
approximately equal. 

Finally it is possible in a strongly unstable case of residual 
oscillations, such as shown in figure 5, that the fuel may 
lose just enough of its energy in splashing to reduce its 
amplitude to where the motion is again smooth. Then, 
because of the instability of the system for small motions, 
the amplitude might again begin to build up. In this way 
continued oscillations of a more or less regular nature would 
occur when the calculations neglecting splashing show un- 
stable motion. This result is important because it shows 
that somewhat irregular small-amplitude oscillations can be 
expected when the ratio of the airplane natural frequency 
to the fuel natural frequency approaches unity, even for 
moderate fuel masses of the order of one-tenth the total 
mass or less. 

EFFECTS OF ASPHERICAL TANKS 

The calculations have been carried out for rigid-body 
motion in spherical tanks only. Actually this assumes that 
for small oscillations the fundamental wave motion in 
spherical tanks approximates rigid-body motion. This 
approximation only applies when the tanks are one-half full 
or less. This restriction is not too serious, however, since; 
the fuel motion will generally have its greatest effect in this 
range. 

It is important to note that the potential energy of the; 
fuels is simply the potential energy of a set of harmonic, 
oscillators located at the positions R f . Thus, the same* 
general analysis will apply whenever the fuel motion in the* 
tank can be represented in terms of harmonic oscillators 
with given effective mass and spring constant. Reference 4 
has already been mentioned as obtaining such a representa- 
tion for the fundamental mode of a rectangular tank 
Usually the fundamental mode will be the most important 
and will involve the greater effective mass. It is conceivable 
that for long tanks the second mode might be of a frequency 
closer to that of the airplane and in that case might be more 
important. In such a case each mode might be represented 
by a separate oscillator. As has been pointed out, the 
damping is mainly due to turbulence and will be more im- 
portant for aspherical tanks. For small motions, however, 
the damping may still be neglected. 

From the general derivation of the equations of motion, 
the most important result was the effect of fuel distribution 
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on the coupling of lateral and longitudinal motions. It is 
plain that these results apply strictly only to spherical tanks. 
Consider for example a tank of triangular plan form located 
in the symmetry plane and oriented symmetrically with 
respect to this plane. Because of the symmetrical orienta- 
tion it can be seen that, although sideslip motion will give 
rise to forward and rearward forces (because the pressure 
forces are normal to the diagonal surfaces), the forward and 
rearward motion will give antisymmetrical lateral forces 
which will cancel. In this case the result would be to feed 
energy from the lateral into the longitudinal motion; this 
condition would be favorable since the longitudinal motion 
is generally well-damped. For an unsymmetrically oriented 
tank of this type, energy could be fed back again from the 
longitudinal to the lateral motion and the problem would be 
quite complicated. 

In general, the results on coupling for spherical tanks 
would be valid to first order for such symmetrical plan forms 
as the rectangular or the diamond-shape ones. For any 
simple symmetrically oriented shape in the plane of sym- 
metry, a loss of energy from the lateral to the longitudinal 
motion might occur. This condition would be favorable. 
Finally, for tanks outside the plane of symmetry the same 
considerations would be valid if the tanks were symmetrically 
placed with respect to the symmetry plane and symmetrically 
shaped with respect to the plane through the tank center 
parallel to the symmetry plane. 


CONCLUSIONS 

The following conclusions may be drawn from the the- 
oretical analysis presented: 

1 . Considerable disturbances of the normal airplane motion 
can be caused by fuel motion. 

(a) The most important factor determining the effect 
is the ratio between fuel and airplane frequencies. When 
these are equal, even moderate amounts of fuel (one- 
tenth the total mass or less) may cause considerable 
disturbances. 

(b) The most usual type of disturbed motion is a some- 
what irregular small-amplitude oscillation and the type 
of motion is strongly dependent on the initial conditions. 

(c) The effects of splashing will be to make the motion 
more stable, and the loss of energy in fuel turbulence may 
make it possible to increase the stability by artificially 
introducing turbulence in the fuel. 

2. The fuel motion may cause coupling between lateral 
and longitudinal motions. 

3. The derivation of the equations of motion for spherical 
tanks may be applied to any tanks where the fuel motion 
may be represented in terms of harmonic oscillators. 

Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 21 , 1950. 
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A STEP-BY-STEP SOLUTION OF THE EQUATIONS OF MOTION IN MATRIX NOTATION 


A matrix method for solving the equations of motion is 
given in reference 9. The first step in this method is to 
reduce equations (19) to a set of first-order equations by 
introducing as new variables Z?</>, D\p, I)£ l} and D{ 2 as de- 
scribed in reference 9 for D(f) and Tty. This transforms the 
equations of motion into a set of nine linear first-order 
equations in the nine variables. In matrix notation the 
equations may then be written as follows when there are no 
applied forces C Y , C n , or C t : 

A(Dq)+Bq= 0 (Al) 

where A and B are ninth-order square matrices and q is the 
column matrix (or vector), the elements (components) of 
which are the nine variables. In partitioned form, 


1= r i 4 j O45 n 

L^54:^4 5 5 J 

r 945; / 4 ”i 

-LcjjseJ 


where / 4 is the identity matrix of fourth order; 0 45 and 0 54 
are zero matrices of order (4X5) and (5X4), respectively; 
and 
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Equation (Al) may be solved for Dq by multiplying 
through by the inverse of A: 

Dq— —A~ l Bq = Pq (A2) 


Equation (A2) shows that, when there are no applied forces, 
the differentiation operator with respect to nondimensional 
time may be replaced by premultiplication with the matrix 
P = —A~ l B. For example, note that 

D 2 q =D(Dq ) =P(Dq) = P 2 q 

and that similar relations would result for higher powers. 

Now, by Taylor’s expansion, the value of q(s-\- As) may 
be obtained from the value of q(s) by the series 

?(«+Aa)= 9 («)+^ Z? 9 (s) + ( -^ )2 D*q(s)+. . . 


Bv use of equation (A2), this equation can be written 

g(s+As)=[/ 9 +^y P+ ( ^ 2 P°-+ . . .] q(8) (A3) 

which is the fundamental recurrence relation used in the 
step-by-step calculation. In equation (A3), 7 9 is the ninth- 
order identity matrix. The set of initial disturbances q 0 
being given, the magnitude of the step As will determine 
the number of powers necessary to obtain a given accuracy 
in the solution. Because of the relatively high fuel fre- 
quencies and because it was desired to obtain the motion to a 
rather large number of periods with reasonable accuracy, 
the series in equation (A2) was used to the sixth power with 
As approximately 1/20 of the airplane period. Thus, the 
matrix relation (A3) was 


where 


q(s+As)=Qq(s) 

= + P + ^f P 2 + • • • 


(A4) 
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